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When inflation is driven by a pseudo-scalar fleld x coupled to vectors as ^ x^F, this coupling may 
lead to a copious production of gauge quanta, which in turns induces non-Gaussian and non-scale 
invariant corrections to curvature perturbations. We point out that this mechanism is generically 
at work in a broad class of inflationary models in supergravity hence providing them with a rich 
set of observational predictions. When the gauge flelds are massless, significant effects on CMB 
scales emerge only for relatively large a. We show that in this regime, the curvature perturbations 
produced at the last stages of inflation have a relatively large amplitude that is of the order of the 
upper bound set by the possible production of primordial black holes by non-Gaussian perturbations. 
On the other hand, within the supergravity framework described in our paper, the gauge fields can 
often acquire a mass through a coupling to additional light scalar fields. Perturbations of these 
fields modulate the duration of inflation, which serves as a source for non-Gaussian perturbations 
of the metric. In this regime, the bounds from primordial black holes are parametrically satisfied 
and non-Gaussianity of the local type can be generated at the observationally interesting level 
/nl~C(10-100). 



I. INTRODUCTION 

In a recent series of papers [T]-[5] a new broad class 
of models of chaotic inflation in supergravity has been 
developed. These models generalize the simplest model 
of this type proposed long ago in ||4j ; see also [St-BU] for 
a partial list of other closely related publications. 

The new class of models [IH3] describes two scalar 
fields, S and <i>, with the superpotential 

W = , (1) 

where /($) is a real holomorphic function such that 
/($) = /($). Any function which can be represented 
by Taylor series with real coefficients has this property. 
The Kahler potential can be chosen to have functional 
form 

K = K{{<^ -^ f,SS). (2) 

In this case, the Kahler potential does not depend on 
(f) = y/2Ke^. Under certain conditions on the Kahler 
potential, inflation occurs along the direction S — Im $ — 
0, and the fleld plays the role of the inflaton fleld with 
the potential 

V[<i>) = \]{<f>l^2)\\ (3) 



This class of models can be further extended [3J |TT] to 
incorporate a KKLT-type construction with strong mo- 
duli stabilization [21], which may have interesting phe- 
nomenological consequences and may provide a simple 
solution of the cosmological moduli and gravitino prob- 
lems pa [23]. 

The generality of the functional form of the inflationary 
potential 1^(</') allows one to describe any combination of 
the parameters and r. Thus, this rather simple class of 
models may describe any set of observational data which 
can be expressed in terms of these two parameters by 
an appropriate choice of the function /($) in the super- 
potential. Meanwhile the choice of the Kahler potential 
controls masses of the flelds orthogonal to the inflatio- 
nary trajectory [TH5]- Reheating in this scenario requires 
considering the scalar- vector coupling ~ cfiF^^F'^'^ [3^, '24] . 
If not only the inflaton but some other scalar fleld has a 
mass much smaller than H during inflation, one may use 
it as a curvaton fleld [25l for generation of non-Gaussian 
perturbations in this class of models [26 ^ 

In this paper, we will study an alternative mechanism 
of production of non-Gaussian perturbations, using an- 
other formulation of this class of models, with 

= i^(($ + $)^5'S'). (4) 

In this case, the Kahler potential does not depend on 
X — \/2Im$, which plays the role of the inflaton fleld 
with the potential 



All scalar flelds have canonical kinetic terms along the 
inflationary trajectory S = Im $ = 0. 



V{x) = \.nx/V2)\'- 



(5) 
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The difference is that now the inflaton field is a pseudo 
scalar, which can have a coupling to vector fields 



(6) 



e^^P" Fpa and a is a diniensionful constant. 



where F^^ 

This coupling is expected to be present since it is compat- 
ible with all the symmetries, including a shift-symmetry 
in X- 

The study of the phenomenological effects of such a 
coupling during inflation has received a lot of attention 
lately [27 -,33j . In particular, it has been shown in [30, 31 
that, if the constant a is large enough, such a coupling 
can lead to a copious production of gauge fields due to the 
time dependence of x- Through their inverse decay into 
inflaton perturbations, these gauge fields yield an addi- 
tional contribution to the scalar power spectrum which 
is both non-Gaussian and violates scale invariance. In 
this way it is possible to obtain non-Gaussian and non- 
scale invariant effects that can be observed by the Planck 
satellite and has not yet been ruled out yet by WMAP, 
although the parameter space corresponding to such a 
signal is relatively small [|33j. In addition, gauge fields 
source tensor modes and lead to a stochastic gravity wave 
signal that could be detected at interferometers such as 
Advanced LIGO or Virgo [32, 34 (see also [35 J. 

Since the new class of infiationary models in supergrav- 
ity needs a coupling between the infiaton and gauge fields 
to have successful reheating, we have to consistently take 
into account the violations of Gaussianity and scale in- 
variance induced by the inverse decay mechanism. This 
is the topic of section |IT] 

A potential threat in this model is the overproduction 
of primordial black holes. As we will see in section |ITTj at 
very small scales, far beyond what is observable by the 
CMB, the produced gauge quanta largely increase the 
curvature power spectrum. At some point, various forms 
of backreaction stops this growth, but by then the power 
spectrum has reached ~ 0(10"^). At such high val- 
ues, a statistical fluctuation might locally increase the 
density so that primordial black holes are formed. In this 
way the non-detection of primordial black holes puts an 
observational upper bound on the power spectrum |36l - 



HT] , which we discuss in section IV Our estimates for the 
late-time power spectrum land a factor of six above this 
bound (compare e.g. (33 1 with (39)). Since we expect our 
estimate to be reliable up to factors of order one, we can- 
not definitively claim that the inverse decay mechanism 
and its interesting phenomenology is incompatible with 
current data, but our result on production of primordial 
black holes highlights a clear tension. 

In section |V] we argue that in supergravity infiation- 
ary models the gauge field can naturally acquire a mass 
through the Higgs mechanism and as consequence of this 



local non-Gaussianity in the observationally interesting 
region '-^ 10 — 100 can be produced via the mecha- 
nism pointed out in [33 . Consider the additional scalar 
field <T with a small mass mg- ^ H during infiation that 
gives mass to the gauge fields via spontaneous symmetry 
breaking. The light field a and, correspondingly, the vec- 
tor field mass, experiences infiationary fiuctuations. In 
the parts of the universe where the value of the vector 
field mass is small, the vector field fluctuations are easily 
produced since the gauge mass quenches the tachyonic 
instability. This in turns leads to a longer stage of in- 
flation because of the additional friction generated by 
the gauge flelds. Meanwhile in the parts of the universe 
where the fluctuations of the light scalar field a make 
this field large, the vector field mass becomes larger and 
infiation is shorter due to the lack of backreaction. As a 
result, fiuctuations of the light scalar field a lead to fiuc- 
tuations of the total number of e-foldings 5N , i.e. to adi- 
abatic perturbations of metric. This effect may generate 
significant primordial local non-Gaussianity without the 
need to require a complicated regime of reheating which 
is necessary in the conventional curvaton scenario. Also, 
in the regime of parameters relevant for this scenario the 
primordial black hole bounds are satisfied parametrically. 

One might suspect that this regime is not very natural, 
as it requires additional fine-tuning of the theory in order 
to introduce the mass with m„ <C H. However, as we are 
going to show, in the context of the supergravity models 
[T]-2] this regime emerges in a very natural way. One 
may consider, for example, the supercurvaton model of 
Ref . [26J , where the role of the curvaton field a is played 
by the radial component of the light complex field S, 
and generalize this model by coupling the field S to the 
vector field. In that case, the vector field acquires the 
mass proportional to a, which experiences infiationary 
fiuctuations. 



In section [VII| we study the evolution of the light field 
a during infiation in our scenario, which is similar to the 
evolution of the curvaton field a in j^HI, so we will con- 
tinue calling this field the curvaton, and use the results 
of for the description of its evolution. In the original 
model of |26| . just as in any other curvaton model [25J, 
adiabatic perturbations of metric are generated by per- 
turbations of the field a after a complicated sequence of 
reheating, expansion of the universe, and the subsequent 
decay of the curvaton field. In our scenario, adiabatic 
perturbations are produced due to the modulation of the 
duration of infiation by the perturbations of the field a. 
As we will demonstrate, this mechanism can easily pro- 
duce local non-Gaussianity in the potentially interesting 
range /nl from 0(10) to O(IOO), even if the coupling 
constant a is not very large. 



Finally, in section |VIII| we find that typical values 
of the coupling constant a considered in this work lead 
to a relatively high perturbative reheating temperature 
T ^ lO^'^GeV. This should be read as a lower limit, since 
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the copious non-perturbative production of gauge fields 
already during inflation could lead to and even higher re- 
heating temperature. This could lead to the cosmological 
gravitino problem [12], but in the class of models with 
strong moduli stabilization and gravitino mass 0(100) 
TeV this problem does not appear |23| . 



II. CMB SCALES: VIOLATIONS OF 
GAUSSIANITY AND SCALE INVARIANCE 



orders in Sx and 6 A one finds 
S = const - / d^xy/^{5x)a {E ■ B) - E ■ B 



(9) 



+ -X6F5F + -SxSFE 



IT ~ 1 

5xSF5F+-{5xfV"'{x) 
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Recently there has been a lot of interest in the effect of 
gauge field production in axion infiation |27H33| . In this 
section we summarize the main points. 

Consider a pseudo-scalar infiaton with a potential suit- 
able for infiation. The symmetries of the theory allow for 
a coupling x^tJ.iyF^'^ to some U{1) gauge sector. This 
coupling is essential for reheating in the supergravity 
models we discussed in section |T1 We will therefore con- 
sider the following bosonic part of the actiorj^ 



S 



lidxf 



:F^ + -XFF + Vix) 



Since all relevant effects arise from the couplings above 
we can safely neglect the gravitational interaction be- 
tween perturbations and work with an unperturbed 
FLRW metric^ We organize the perturbation theory 
based on the equations that we are able to solve. Con- 
sider two classicaj^fields A{x, t) and that solve these 
two coupled differential equations 

X + 2,Hx + V'{x)^(^{E-B), (7) 
V^l-ax'V X 1=0, (8) 

where E = -A/a, B = a^^y x A tmd E ■ B = -FF/A 
are computed from A. 

Now let us look at the action expanded around x and 
A, i.e. S[x + Sx,A + SA]. Organizing the result at various 



Notice that in the existing hterature, such a coupUng is usually 
associated with interaction of the axion field with vector fields, 
with a coupling In our approach it is not necessary to 

associate the pseudo scalar field with the axion field with the 
radius of the potential ~ /, so we normalize the coupling in 
terms of the reduced Planck mass Mp, which we then set to one, 
and consider the following interaction term — ^X^m"^''''- 
We are neglecting vector and tensor modes and the slow-roll sup- 
pressed interactions coming from the solution of the constraints 
on the lapse and the shift. 

Here we are assuming that the occupation number of the relevant 
gauge modes is large enough that one can approximate the result- 
ing electromagnetic field with a classical one. This assumption 
is implicit in all other approaches so far. 



where again the classical background fields x and A solve 
([t]) and (|8|. Notice that there is a "tadpole" for Sx due 
to the fact that at the background level we solved an 
inhomogeneous equation for A but just a homogeneous 
one for x- From this term one also sees that Sx will source 
SA'^, hence it will modify the constraint. The equations 
of motion in Coulomb gauge diA'' = are 



aSAj 



dliSA.) 



aHSAi - axV x (SA) 



(10) 



aSx"^ X A-a{WSx) x 1- dt{ad,{SA'^)) , 
iSx)+iHSx~V^Sx + V"Sx^ (11) 
^(^{FF)-FF-2Ff) , 

ad^d^iSAf = -aW{Sx) • V x 1. (12) 



The solution for the constraint equation for SA'^ is 



SA"{x,t) = a-^ / d^y 



^3 aV{5x)-VxA 



A'k\x - 



V\ 



(13) 



Unfortunately this coupled system of equations is hard to 
solve. Hence [30', "ST] made the approximation of neglect- 
ing all terms quadratic or higher in Sx, SA and A. This 
is a good approximation as long as EE (or equivalently 
{E ■ B)) is not too large (a more quantitative condition 
is given in ( |29[ )), which is the regime we will discuss in 
this section. In the next section we will see that, since 
{E ■ B) grows with time, towards the end of infiation this 
description in not valid anymore, and one has to take 
backreaction into account. 



Solving the approximated equations of motion 



aAi - 



(14) 



aUAi - axbV x A = 
Sx + "iUSx - y^Sx + V"Sx ^a{{E ■ B) - E ■ b)(15) 

one finds a tachyonic enhancement of the gauge fields. 
For the growing mode of one of the two polarizations of 
the gauge field we get 



A 



1 



e</2w^_,.i/2(2zfcT). 



(16) 
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Here W\,ii(x) denotes the Whittaker function, and ^ is 
defined aall 



'2H' 



(17) 



As we see, the relation between the coupHng constant a 
and the value of ^ 60 e-foldings before the end of infla- 
tion is model dependent, but for our model there is an 
approximate relation which is valid for the parameters 
that we are going to explore: 



a ~ 15^. 



(18) 



For ^ > 1 the new coupling therefore leads to gener- 
ation of perturbations of the vector fields around hori- 
zon scales. The produced gauge fields then change the 
dynamics of x ^^nd H. The cosmological homogeneous 
Klein-Gordon equation and the Friedmann equation get 
extra contributions from the gauge fields and can now be 
written as 



3Hx 
1 



dV 

dx 

2 . V 



a{E- B 
1 



They are computed as 

{E-B)^-\- 
,E^ + B\ 1 



dke^\A\^ 

dr 



47r2a4 



dkk'^ 



(19) 
(20) 

(21) 
(22) 



After renormalization, one can reduce the integration in- 
terval to the region < < 2^, which is where the 
enhancement in the (derivative of the) gauge field takes 
place. 



From the homogeneous Klein-Gordon equation ( 15 ) 



one reads off that the infiuence of the produced gauge 
fields on the homogeneous dynamics of x ^^nd H can be 
safely neglected as long as 



a{E-B) 
SHx 



« 1, 



3IP 



< 1. 



(23) 



* Note that we have some minus signs different from 1301 . but this 
is a matter of conventions. We will work with a model that has 
X < during inflation and deflne ^ to be positive. The sign of 
{E ■ B) is always opposite to the sign of x- Therefore the physical 
effect of the tachyonic enhancement is always that inflation is 
prolonged. To be precise: when x is negative, the growing field 
is actually the opposite polarization, i.e. A_, which makes that 
{E ■ B) > (see, for example, equation (8) in [28]). 



Of these two conditions the first one is always the most 
stringent. When it stops to hold, backreaction on the 
homogeneous evolution becomes important and the evo- 
lution of X E^nd H will be slowed down, which makes in- 
fiation lasts longer. We will see in the next section that 
backreaction on the inhomogeneous equation for 5x hap- 
pens even earlier. In this section we focus on the regime 
in which all of these effects are negligible, which e.g. for 
a quadratic potential corresponds roughly to ^ < 4. This 
is appropriate for the description of CMB scales. 

Now we move to the power spectrum. The copiously 
generated gauge fields may, by inverse decay, produce 
additional perturbations of the infiaton field 5x, propor- 
tional to the square of the vector field perturbations. As 
was shown in [3(3131], this can be described (up to backre- 
action effects to be described in the next section) by using 
(151. The inclusion of the source term leads to an extra 



contribution to the power spectrum of the curvature per- 
turbation on uniform density hypersurfaces Q = —^5x, 
which has been computed in |30|, I31| (we present a quick 
estimate in Appendix ([B])) 

A2(fc) ^ ^i Jk) (1 + ,^(fc) hiOe'^^) , (24) 

where /2 {£) was defined in |30l |SI] and can be computed 
numerically (a useful large ^ approximation is given in 
(|bT4|) and 



27r|x| 



(25) 



is the amplitude of the vacuum infiationary perturbations 
as in standard slow-roll infiaton. WMAP [43 has mea- 
sured A^ siX*^*) = 2.43 • 10-9, where = 0.'002Mpc"^ 
is the pivot scale that we will take to correspond with 
iV = 60 e-foldings before the end of infiation. The second 



term in ( 24 1 violates both scale invariance (and Gaussian- 
ity as we will see below), since it comes schematically 
from [SAY, i.e. the square of a Gaussian which grows 



with time as in (16 1 



We move to the bispectrum. The produced gauge fields 
lead to equilateral non-Gaussianity in the CMB |30l EI] 



Tnl 



/3(e), 



(26) 



where /a (^ ) another function defined in |30[ [3T] , which 
can be computed numerically (see (D7l for a useful ap- 



proximation). The amount of non-Gaussianity, therefore, 
depends exponentially on S^. Between ^ = and ^ = 3 it 
grows from 0(1) to 0(10"*) and most of the growth takes 
place in a small interval around ^ ~ 2.5. 

The analysis of [25] showed that the bounds coming 
from the power spectrum (especially from WMAP plus 



5 



14 
12 
10 
8 
6 
4 
2 



10 



20 



30 



40 



50 



60 



N 



FIG. 1. The evolution of the infiaton field x, as a function of the 
number of e-folds N left to the end of inflation (time is moving 
to the left) for g[7V = 60] = 2.2. The result in dashed blue does 
take backreaction from the sources in equations l |19| and | |20[ into 
account, the result in red does not. It is clear that backreaction 
prolongs inflation. 




FIG. 2. The evolution of the Hubble scale as a function of A' for 
^[N = 60] = 2.2. Again the dashed blue line is the result corrected 
for backreaction from the sources in equations \19\ and l |20[ l. 



ACT, because of the violation of scale invariance) and 
from the bispectrum (from WIVIAP) are compatible, with 
the former being typically slightly more stringent. Spe- 
cifying a confidence region in ^ requires assuming some 
prior for this parameter. The physically best motivated 
prior is log-flat in ^ reflecting the fact the scale of the 
dimension flve coupling x^P could be anywhere (with 
strong indications that it should be below the Planck 
scale [44J). In this case at 95% CL one finds ^ < 2.2. A 
flat prior on ^ leads to ^ < 2.4. 



III. VERY SMALL SCALES: STRONG 
BACKREACTION 



In this section we want to estimate the power spectrum 
and bispectrum towards the end of inflation, i.e. on scales 
that are too small to be observed in the CIVIB. The only 
observational handle available in this regime is the non- 
detection of primordial black holes, which puts an upper 



bound on the power spectrum 



To make these estimates it is essential to recognize that 
many of the formulae described in the previous section 
and given in the literature about inverse decay are valid 
only in the regime in which backreaction on the inhomo- 



geneous equation for dx is small (see (29l). As we show 



in the following, the scales relevant for the production 
of primordial black holes leave the horizon when backre- 
action is large. The authors of [45j did not account for 
backreaction and therefore their conclusion that gauge 
field production during inflation leads to black hole pro- 
duction might be premature. 

For concreteness, we will consider a quadratic potential 
^(x) — ^Tn^X^i with the mass chosen such that at the 
pivot scale fc^ (that we take to correspond with N = 60) 
we get A2(fc^) = 2.43 • 10"^ 

Let us first look at the dynamics of x and H. As we 
already discussed, when enough gauge field quanta have 
been produced, the conditions in ( p3| stop to hold (the 
inequality for {E ■ B) is violated first) and x and H are 
slowed down. As a result, inflation lasts longer. Let us 
check this. The behavior of x, H and as functions of N 
(the number of e-folds left to the end of inflation) follows 
from simultaneously solving (17), ( [T9| and (20). In Fig- 
ures ([T]) and ([2| we have plotted the solutions for x(A^) 
and hIn) , with and without backreaction taken into ac- 
count. For £^{N = 60) = 2.2, the effect of backreaction 
becomes 10% around iV = 11. 

Now let us consider perturbations. Of course they will 
be affected by the backreaction on the homogeneous dy- 
namics X and H that we described above, but there is 
more. Let us consider (10l-(12l. In the last section we 



solved for A in a homogeneous background and used that 
result (16) to compute the source term in the equation 
for X perturbations. But as SA and dx grow larger to- 
ward the end of inflation (both of them grow as e^'^^) 



the source in the right-hand side of (10) can not be ne 



glected anymore. If we were able to solve this equation, 
we would find that E-B now depends on the perturbation 
dx- Expanding E ■ B, which is the source term in (111, 



in powers of 6x we would find several new terms includ- 
ing additional friction and a modified speed of sound. In 
|28l [35] it was proposed how to estimate these effects in 
the regime of strong backreaction by just considering the 
additional friction term 6x- The equation of motion for 
the perturbation dx becomes 



6x + iPHSx 



-Sx 



Sx ■ 



E ■ B -{E- B) 



with the additional friction term 



1 - 27rCa- 



(E-B) 
3Hx 



(27) 



(28) 
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Here the new term in /3 is caused by the dependence of 
{E ■ B) on X (via its dependence on ^) . The behavior of 
/3 has been plotted in Figure ([3|. It is always positiv^ 

The new source of backreaction can be neglected as 
long as 



< 1. 



(29) 



Note (from comparison with (23)) that the factor of 27r^ 
makes that backreaction on the power spectrum will be- 
come significant before backreaction on H and x does. 
For £^{N = 60) = 2.2 we find that backreaction becomes 
of order 10% [p = 1.1) at TV = 22. 



The modified equation of motion (27 1 suggests that (as 
was already noted in |32|i see also appendix (|B|) we can 
estimate 



5x 



a{E-B-{E-B)^ 



which leads to the power spectrum 



(30) 



A^(fc) ^ (C(x)2) 



^{E ■ B) 

'WW 



(31) 



This estimate turns out to be particularly good in the 
regime in which we can check it, i.e. when ^ ^ 4 when 
the backreaction is negligible and we can compare with 
(24 1 ((see appendix (|b]))). This gives us confidence to 
use it also in the strong backreaction regime. It is easy 
to see that when backreaction becomes large, the second 



term in (28 I dominates, and we end up with 



p-1 
100^ 
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FIG. 3. Evolution of 1) as function of N, for i{N = 60) = 2.2. 
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FIG. 4. Evolution of the power s pect rum as function of A^, for 
^(N = 60) = 2.2. The expression |24]l that does not take back- 
reaction into account is in tinily dashed blue. In solid red is the 
estimate l |31| . When backreaction becomes significant this esti- 
mate coincides with the late-time estimate (27r^[A''])~^, in largely 
dashed green. 



(32) 



IV. BOUNDS FROM PRIMORDIAL BLACK 
HOLES 



The estimate ( |31| for the power spectrum has been 
plotted in Figure (|4| together with the formula (24 1, 
valid only when backreaction is negligible. Indeed, in 
the regime of strong backreaction the power spectrum 
asymptotes the estimate in ( 32 1 . At the end of inflation 
we have f 



6.7 (for ^{N = 60) = 2.2), which gives 



7.5 • 10" 



(33) 



^ We work with negative x which yields positive {E ■ B), while 
working with X > gives {E ■ B) < 0. 



Now let us try to compare this with the existing bounds 
on the power spectrum coming from the non-detection 
of primordial black holes. These will form if at hori- 
zon re-entry (i.e. smoothing C on scales of order H) we 
have C > Cc, with Cc ^ 1 denoting the critical value 
leading to black hole formation. If one assumes that ( 
follows a Gaussian distribution (with (C) = 0) one can 
express the probability of having C > in terms of the 
variance (('■^) by analyzing the Gaussian probability dis- 
tribution function. This probability corresponds to the 
fraction of space /3 that can collapse to form horizon-sized 
black holes. Hawking evaporation and present day grav- 
itational effects constrain this fraction /3. Typically one 
finds (3 in the range (lO^^^ — 10~^), with the strongest 
bounds coming from CMB anisotropics |37| (spectral dis- 
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tortion and photodissociation of deuterium lead to a 
bound /3 < 10~^°, as for example in [M]). Setting Cc = 1 
gives for the upper bound on the power spectrum [JD] 



A^^(fc) ~ (C(a;)^) ~ 0.008 - 0.05. 



(34) 



Here the lower bound corresponds to /3 = 10 and the 
upper bound to /3 = 10^^. 

However, in our case C does not follow a Gaussian dis- 
tribution. Instead we have (see Appendix ([B|) 



c = - 



a(^E-B- {E-B] 



(35) 



The stochastic properties of the vector field 5 A are close 
to those in a free theory, i.e. it has Gaussian perturba- 
tions around {5A)—Q. As a consequence we can write 



El 



C = .9^ - (.9^) 



(36) 



with g a Gaussian distributed field. This model was stud- 
ied in HO] and we follow that derivation (see also |551[55] 1. 
The probability distribution function of C, follows from 
setting P{()d( — P{g)dg, and takes the form 



PiC) 



V27r(C a2)a 



(37) 



with = (g^). For a given value of /3 we can again 
infer a^. Setting t = -i^ + 1 (and tc = ^ + 1) we have 
d( = a^dt which gives 



At the end of inflation, the total mass concentrated 
in the volume associated with perturbations leaving the 
horizon N e-foldings before the end of inflation with the 
Hubble constant H can be estimated by 



M 



N 



4 3 4^Af2 



(40) 



where we reinserted the reduced Planck mass Mp, which 
was set to one in the rest of the paper, and H is cal- 
culated at the end of inflation. In order to study the 
subsequent evolution of mater in the comoving volume 
corresponding to this part of the universe, one should 
distinguish between two specific possibilities depending 
on the dynamics of reheating after infiation, discussed in 
section IVIIII 

If reheating is not very efficient, then the universe for a 
long time remains dominated by scalar field oscillations, 
with the average equation of state p = 0. In this case, 
the total mass in the comoving volume does not change, 
and therefore at the moment when the black hole forms, 
its mass Mbh is equal to Mm evaluated in (40 1. For 



the parameters of our model, this gives an estimate (see 
appendix [e| for details) 



M 



BH 



10 e 



(41) 



On the other hand, if reheating is efficient, then the 
post-inflationary universe is populated by ultra relativis- 
tic particles and the energy density in comoving volume 
scales inversely proportional to the expansion of the uni- 
verse. In this case, the black hole mass can be estimated 
as (see appendix IeI) 



Msff ~ lOe^^ g 



(42) 



/ P(C)dC 

'Cc 



-.dt =Erfc 



(V¥) ' (38) 



where Erfc(a;) = 1 — Erf(a::) is the complementary error 
function. Taking again j3 in the range 10~^^ — 10~^ one 
gets a tighter upper bound on the power spectrum than 
in the Gaussian case: 



~ 1.3 • 10""^ -5.8- 10"^ 



(39) 



Now let us estimate what value of (5 is relevant for our 
investigation. 



In our estimates of the black hole production we will 
assume the latter possibility, though in general one may 
have a sequence of the flrst and the second regime. The 
flnal conclusion will only mildly depend on the choice 
between these two possibilities. 

Now, the bounds on /3 in terms of the would-be black 
hole mass Mbh were given in |36l and updated in fHTj . 
Here we follow the result in fSTJ F] Using ( 38 1 and our 
estimates of the black hole mass as a function of N , we 
can translate this into a bound on the power spectrum 
as a function of N . The result is in Figure [5] 



Our estimate (33 1 violates this hound for all A'^ < 20 by 



a factor of about six. Since we have made some approx- 
imations both in deriving the late time power spectrum 



^ Here we can safely neglect the linear term, which is just the stan- 
dard vacuum slow-roll contribution to C,. See also our estimate 
for fiqi^ at small scales in Appendix |[D]|. 



^ However, we do not take the constraints for Mgu < 10** g into 
account, as these are either very model dependent, or assume 
that black hole evaporation leaves stable relics. 
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N 



the expansion history of the universe and therefore speci- 
fies NcMB — 64, for the N corresponding to CMB scales 
(see appendix [E| for a derivation). This is required for 
consistency but changes the numerics very Httle. There- 
fore in all other sections we still use Ncmb — 60. 

For the matter domination regime, the black hole 
masses would be greater, for a given iV, see (41 1, and 



therefore we would have a slightly stronger constraint on 
^ and a. We find C ^ 1-3 which corresponds to a < 20. 
Instead of concentrating on it, we will now investigate the 
model where non-Gaussian perturbations may be gener- 
ated for much smaller ^ and a, without leading to the 
primordial black hole problem. 



FIG. 5. Evolution of our estimate for the power spectrum as a 
function of A'^. In dashed red is the result for ^[A'^ = 64] = 2.2. 
Other lines are for g[7V = 64] = 2.5 (solid brown), ^[AT = 64] = 2 
(solid blue), £_[N = 64] = 1.5 (solid green), i[N = 64] = 1 (solid 
yellow) and ^[A'^ = 64] = 0.5 (solid orange). The black hole bound 
is in dashed black. 



LOCAL NON-GAUSSIANITY FROM HEAVY 
VECTOR FIELDS 



and in deriving its observational upper bound, our es- 
timate could well be off by some order one factor and 
therefore we can not draw a definitive conclusion. It is 
clear though that the parameter values giving rise to an 
observable but not yet ruled out violation of scale invari- 
ance and non-Gaussianity in the CMB-window produce 
a late time power spectrum that comes at least very close 
to the primordial black hole bound. A more precise com- 
putation is needed to establish whether this bound is 
actually violated or not. 

However, if such a computation revealed that primor- 
dial black holes do indeed constrain these models, that 
would yield a much stronger bound on ^ as the ones com- 
ing from non-Gaussianity and the violation of scale in- 
variance. Since we have seen that the power spectrum 
has a late-time asymptotic of (27r^[A^])~^, this problem 
persists on a wide range of values for ^. 

For all values of ^, our estimate for the power spectrum 
sharply increases before the end of inflation, the closer to 
the end the smaller <^ is. However, if we disregard black 
hole bounds for M^h S 10^ gj which rely on uncertain 
model dependent assumptions, there are no black hole 
bounds for < 8. From figure [s] we then see that we get 



CMB. 



< 1.5 



(43) 



for the bound on f at CMB scales from primordial black 
hole production. In terms of the coupling constant a, 
this bound implies the constraint 



5, 23. 



(44) 



This bound is derived using (42 1, i.e. radiation domina- 



Now let us turn to a scenario, described in [331, in 
which the produced gauge fields are massive. The pro- 
duction of gauge quanta decreases with the mass of the 
gauge fields: for tua ~ £,H all production is killed. In this 
scenario, the gauge fields get their mass via the Higgs 
mechanism. Fluctuations in the Higgs field h lead to 
fiuctuations in niA, which in turn generate fluctuations 
in the amount of produced gauge quanta, and therefore 
in the amount of extra friction in the dynamics of x E^nd 
H. In the end, one has perturbations in AiV, namely 
the number of extra e-folds of inflation due to gauge field 
production. This leads to a non-Gaussian signal in the 
CMB of the local type [33J . Using the 6N formalism one 
finds 



/■local 



102 



Cio-3 



(45) 



tion right after the end of infiation. This assumption fixes 



Here h is the Higgs-like field responsible for the sponta- 
neous symmetry breaking that gives a mass to the gauge 
fields, e is its U{1) charge, itia = eh and assumed a 
quadratic potential H — 

For a complete description we refer the reader to the 
original reference [33', section 7. Here we only want to 
stress that this scenario can also work for ^ ~ 1. Then 
it will surely satisfy the bounds from primordial black 
holes. However, this scenario requires that m/i <C H, 
which seems rather unnatural. In the next two sections 
we want to explain how this model, that can produce an 
observable local non-Gaussian signal, can be embedded 
in a supergravity model that is a slight extension of the 
models considered so far. In this model one can easily 
have m/i <C as we will argue. 
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VI. CURVATON SCENARIO AND CHAOTIC 
INFLATION IN SUPERGRAVITY 



We go back to the supergravity scenario described in 
section [T] that we want to extend by making the vector 
fields heavy because of the Higgs effect. We will start 
with the simplest version of chaotic inflation proposed in 
[4], with the superpotential 



and Kahler potential 



(46) 



(47) 



The Kahler potential does not depend on the phase of 
the field S and on the imaginary part of the field $. 
It is convenient to represent the fields S and ^ as S — 
a e'^/V2 and <f> = {(j) + ix)/V2. The field x plays the role 
of the inflaton field, with the quadratic potential, as in 
the simplest version of the chaotic inflation scenario : 



Note that = for xcr <C 1. During inflation 
rri^ <C H^, and therefore inflationary perturbations of 
the fleld a can be generated. At xcr > 1, the effec- 
tive mass squared of the field a is dominated by the 
term |m^x^f^ — ^H^cr'^ > m'^ ■ For cr ^ 1 one still 
has ^ H^, so the perturbations of the field a are 
generated in this regime as well. At cr > 1 the poten- 
tial becomes exponentially steep, ^ H^. Therefore 
infiationary fluctuations of this field are generated only 
for a < 1. The steepness of the curvaton potential at 
CT > 1 protects us from extremely large perturbations 
of the curvaton field which otherwise could be produced 
during eternal infiation |47l H5] . 

If one does not take into account the curvaton fiuctu- 
ations in this scenario and study only the usual infiaton 
fluctuations J^\, then the COBE normalization requires 
m ^ 6.35 X 10~^, in the system of units Mp = 1. With 
an account taken of vector fleld production modulated 
by curvaton fluctuations, inflation lasts a bit longer, and 
therefore one has to take the inflaton mass a bit larger. 
For ^ = 1 one needs m ^ 6.45 x 10^^. 



v{x) 



-X 



(48) 



Near the inflationary trajectory with 5 = 0, the mass 
squared of the real part of the fleld <i> is to| = QH^ + 
rn? , where — rv?x^ j^S. During inflation to| > 6i/^. 
Therefore the fleld quickly rolls towards its equilibrium 
position at = 0, and perturbations of this fleld are not 
generated. Meanwhile the fleld 5 is light, its mass is well 
below _ff, and can play the role of the Higgs fleld discussed 
in the previous section. Hence we assume that 5 has 
some U(X) charge e and a minimal coupling to gauge 
flelds 1(9^ + ieA^)5p. This is the main new ingredient 
that we introduce now to our model. 

Both components of the fleld 5 may remain light dur- 
ing inflation, and therefore inflationary perturbations of 
these flelds can be generated [TT] . The fluctuations of the 
fleld a are generated just as the fluctuations of the curva- 
ton fleld [25, .26j , though in the present context its fluc- 
tuations will be important for a different reason. They 
will be responsible for the fluctuations of the mass of the 
vector fleld, in accordance with the scenario described in 
the previous section. 

The potential of the flelds x, cr at = is 



VII. STOCHASTIC APPROACH 



We will begin our study with investigation of the be- 
havior of the distribution of the fluctuations of the cur- 
vaton fleld (T, following |2^. During inflation, the long- 
wavelength distribution of this fleld generated at the 
early stages of inflation behaves as a nearly homogeneous 
classical fleld, which satisfles the equation 



ma + K = 
which can be also written as 



(52) 



(53) 



However, each time interval new fluctuations of the 
scalar fleld are generated, with an average amplitude 
squarecF] 



Ml 

2^ ■ 



(54) 



For (T <C 1 one has 



(49) 



9 9 ? ? 9 9 4 

y{x.o) = — + — + • (50) 



The effective mass squared of the fleld a in this regime is 



2 2 2 2 2 

m_ = m H — m y a 

a 4 



(51) 



* For a real massless field we would get (<5a"^) = ■ An extra 
coefliicient 2 appears in | |54| because the field S is complex, so its 
absolute value changes faster because of independent fluctuations 
of its two components. One could argue that in the unitary gauge 
we only have one scalar degree of freedom. However, unitary 
gauge is problematic in the description of the Brownian motion 
and cosmic string formation in the Higgs model. We present 
the results which should be valid in the regime of small gauge 
coupling constant e. Our main conclusions are unaffected by this 
factor of 2 issue. 
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The wavelength of these fluctuations is rapidly 
stretched by inflation. This effect increases the average 
value of the squared of the classical field cr in a process 
similar to the Brownian motion. As a result, the square 
of the field a at any given point with an account taken 
of inflationary fluctuations changes, in average, with the 
speed which differs from the predictions of the classical 
equation of motion by : 



dt 



2Vaa 
3H 



27r2 



(55) 



Using 3i?x = ~^x' o'^^ ^^'^ rewrite this equation as 



(56) 



By solving this equation with different boundary condi- 
tions, one can find the most probable value of the locally 
homogeneous field a. In the supergravity model that we 
study, this equation has an important advantage |26| : 
For a very broad choice of initial conditions, its solutions 
rapidly converge to an attractor solution 



j^2/3 ^ 



• r(i/3) \ 
2\/3^2/3 Vr(2/3)y 



1/2 



0.19 m^/^^. (57) 



If the amplitude of perturbations of metric is determined 
mainly by the infiaton perturbations, one should have m 
in the range of 6.35 x 10~^ to 7.2 x 10~^, depending on 
^ (these two values correspond to ^ = and ^[A^ = 60] — 
2.2, respectively.) 



Now we will use the expression for /nl*^' (|45| obtained 



mA 



33], for mA 



ea and H 



mx 
76' 



- 2.15 X 10 
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under condition 



(58) 



As we see, this expression is very sensitive to the value of 
the gauge coupling e. For any value of ^ and m, one can 
calculate AAmax, and then one can find a particular value 
of the gauge coupling e which results in the perturbations 
with a desirable value of /^l^': 



Any further decrease of f well below 0.5 rapidly shuts 
down generation of the vector field fiuctuations and 
makes AA^max exponentially small, which does not allow 
production of large non-Gaussianity. One could hope, 
following ( 58 1 , that it is possible to compensate any de- 



crease of ^ by an increase of e, but one cannot increase 
e too much without violating the condition mA ^ 
required for the vector field production. 

Until now, we assumed that the amplitude of the met- 
ric perturbations is dominated by the usual infiaton per- 
turbations. However, one may consider models with 
TO ^ 6 X 10^^ and study a possibility that the leading 
contribution to perturbations of metric is given by the 
curvaton-related perturbations. As an example, we will 
assume that m ~ 10^^, in which case infiaton perturba- 
tions give a negligibly small contribution to the COBE 
normalized amplitude of perturbations of metric. Us- 
ing (7.23) of [33] and (57 1 one finds that in this case 
COBE normalization for ^ = 1 requires e w 8 x lO""*, 
which yields ~ 7, which is almost too small to be 
observed. By lowering ^, one can exponentially decrease 
AA'iiiax and get non-Gaussianity in the observable range 
for somewhat larger values of the gauge coupling. 

Now let us add a bit of black magic. When we calcu- 



lated o'(x) in (57 1, in fact we calculated the typical devi- 
ation of the field a from 0, or, in other words, the width 
of its distribution. The field a can be twice as large as 
our result indicates, even though the probability of that 
would be somewhat suppressed. More importantly, it can 
be much smaller than the value which we calculated; the 
probability of this to happen is nearly fiatly distributed 
in the interval from to cr{x)- 

Notice that the amplitude of perturbations of the met- 
ric is proportional to e^a, which we should keep constant 
to maintain the COBE normalization. Meanwhile /nl is 
proportional to e^cr^ = Thus, for the given am- 

plitude of perturbations of metric, /nl ^ l/f, so one 
can get /nl which is two times greater than our estimate 
in those parts of the universe where a = cr(x)/2. As we 
just argued, parts of the universe with a < a{x) are quite 
abundant. Thus we can easily have a proper amplitude 
of perturbations and /nl in the range of 0(10) or higher 
in this scenario. 



e ~ 6 X 10" 



1/9 



/ j'local\ 

^^Unl ) 



1/6 



(59) 



One may consider some numerical examples. For ^ ~ 1 
one finds AA'max ~ 1 (the example considered in [331). 
Then for m ~ 6.5 x 10~^ one may have, for example, 
f^l""^ - 30 for e - 2.8 X lO^^. For ^ - 0.5 one finds 
AA^max ^ 0.044, so one may have /^°£^^ ~ 30 for e ~ 
6.7 X 10~^. Thus the required gauge couplings are pretty 
small. They could correspond, e.g., to gauge interactions 
with vector fields in the hidden sector. 



One can make a similar trick in the previously studied 
regime where the amplitude of perturbations of metric 
was dominated by the infiaton contribution. In that case, 
this amplitude does not depend on cr, whereas /nl is 
proportional to a^. This means that one can live, with 
a significant probability, in the where /nl is much 

smaller than the result obtained above. 

Thus, just like in the general curvaton scenario, we 
are dealing with the "curvaton web", where the universe 
becomes divided into different exponentially large parts 
with different values of the amplitude of perturbations of 
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metric and different values of /nl ITT] . 



VIII. GAUGE FIELD PRODUCTION IN SUGRA 
INFLATION: REHEATING 



We have found that the couphng needed for re- 

heating in (the pseudo scalar variant of) the new class 
of SUGRA inflation models proposed in [1-3 can as well 
yield an observable non-Gaussian signal. It only remains 
to be seen what the effects of the typically needed values 
for ^ are for the reheating in the combined model. 

In [3] the reheating temperature Tji for the decay of a 
scalar inflaton field to two photons due to the interaction 
j(I)F'^ was estimated as 

Tr w V2a X lO^GeV. (60) 

A similar estimate is valid in our case. One may also 
represent it in an equivalent way using relation j = — , 

. 2 

and an expression for the slow-roll parameter e — ^rp- 



TRKi -^x lO^GeV. (61) 



As long as one can describe reheating particle by 
particle decay, reheating in inflationary models of this 
type does not depend much on whether the inflaton field 
is a scalar or a pseudo scalar. In both types of models, 
one may consider interactions with a <C 1, which results 
in reheating temperature Tr < 10^ GeV. This solves the 
cosmological gravitino problem for gravitino in the typi- 
cal mass range 7713/2 ^ 1 TeV. 

However, for a > 1, which is required for production of 
non-Gaussianity in the models based on the pseudo scalar 
inflaton, an estimate described above gives Tr > 10^ 
GeV. It is good for the theory of leptogenesis, but it could 
be bad from the point of view of the gravitino problem. 
Moreover, for a > 1 an entirely different mechanism of 
reheating is operating. At the end of inflation, when the 
time-dependent parameter ^ grows and becomes large, 
a significant fraction of the energy of the infiaton field 
gradually becomes converted to the energy of the vector 
field (see Figure [6]) . This is a very efficient mechanism, 
which may lead to a very rapid thermalization of energy 
in the hidden sector. This may exacerbate the gravitino 
problem. Fortunately, this problem does not appear for 
superheavy gravitino with mass 703/2 ^ 10^ TeV. Such 
gravitinos appear in many versions of the models of mini- 
split supersymmetry, which became quite popular during 
the last few years, see [50J and references therein. 



1^ 
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FIG. 6. Evolution of the normalized energy of the vector field, x = 
I (S2 + B^) /3H'^ as a function of N, for ^[N = 60] = 2.2 (solid 
red), ^[AT = 60] = 1.0 (largely dashed blue) and ^[Ar = 60] = 0.5 
(tinily dashed green) . 



IX. CONCLUSIONS 



The new class of chaotic inflation models in supergrav- 
ity needs a gauge-gauge-inflaton coupling for reheating. 
The inclusion of this coupling can produce gauge fields 
and can provide a Planck-observable but not yet ruled 
out non-Gaussian signal in the CMB. 

In this article we have studied two possible realizations 
of this scenario. Taking the parameter ^ ~ 2.2 — 2.5 {a ~ 
32—37) produces a large amount of gauge quanta, that by 
inverse decay give rise to an equilateral non-Gaussianity 
in the CMB, as studied in |30t I31|. However, we have 
estimated that towards the end of inflation the power 
spectrum grows so much that the model may be ruled out 
because it overproduces primordial black holes. As our 
order-one estimate lands within a factor of six from the 
critical black hole bound on the power spectrum (with the 
non-Gaussian nature of the signal taken into account), 
we need a more precise computation to draw a definitive 
conclusion. 

In the second scenario, where the produced gauge fields 
are massive due to the Higgs effect in presence of a light 
curvaton-type field, one can take a smaller value for ^, 
of order 0.5 - 1, corresponding to a from 8 tol5. Then 
the model is free of black hole trouble. In this case, fluc- 
tuations in the curvaton field modulate the duration of 
infiation and can give rise to adiabatic non-Gaussian per- 
turbations of the local type with /nl ~ 0(10 — 100). For 
smaller values of a, we return to the standard chaotic in- 
fiation scenario with Gaussian adiabatic perturbations. 
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In j30| 131] the power spectrum ( 24 1 has been obtained 
by the Green's function method. In[32, a quick estimate 
was introduced to compute the power spectrum in the 
case of large backreaction (/3 ^ 1). Here we want to re- 
view and further explore this estimate, showing how it 
leads to (31 1 and also how, in the case of negligible back- 



reaction, it approximates the precise result (24 1 within a 
factor of two. 



Appendix A: Variance of E ■ B 



The full equation of motion for the perturbation Sx is 
(in real space) 



The variance of £^ • _B is defined as 

= {{E-Bf) - {E-Bf 
^{E,E,){B,B,) + {E.,B,){B,E,) 



(Al) 
(A2) 



Sx + WH5x 



with 



-5x + rn?5x — ct 



E- B-IE -B) 



(Bl) 



We find 



{E.,E,){B,B,) 



{E,B,){B,E,) = 



dkdq 
(2^ 



A'{k)\'\Aiq)\'q*k' 



(fflk(fD,qei{k)ei{q)e*{k)e*{q) , 



dkdq 
(2^ 



A{k)A'*{k)A'{q)A*{q) q^k^ 



d^^lkd^^lqe,{k)e,{q)e*{k)e*{q). 



(A3) 



The angular integral gives (47r)^/3, i.e. a third of the 
whole sphere. The integrals over the modulus are similar 
to the one in [31 and are computed in the same way 



^/^K(mV.2.2.10 



e— ,(A4) 



1 



dk 



a4 J (27r)' 



-\A{k)\''k^ ~ 1.9 • 10-5^e2"« . (A6) 



Putting things together one finds 



(47r)2 



F74 

4 „2TrC 



-2.0-10 -^e 



-{E-B) 



(A7) 
(A8) 



1 - 27r^a- 



(E-B) 
iHx 



(B2) 



Near horizon crossing we can estimate d ^ H. Since we 
have, near horizon crossing, — the first term can- 
cels the third one. The second term can be approximated 
as 3f3H^Sx- The last term on the left hand side is just a 
slow-roll correction and can be discarded. This directly 



gives 



Sx 



[e-B - {E-B) 



3/3i/2 



(B3) 



and therefore we have 



H 



C^-^Sx^ 
X 



a(^E ■ B - {E ■ B)^ 



il3Hx 



(B4) 



For the position space two point function of ^ we imme- 
diately get 



(C(x) 



ff2 



<5x^> 



( a{E-B} y 
\ 3f3Hx j 



(B5) 



with a the variance computed in the previous subsection. 
To compare the position space power spectrum with 



13 



the momentum space power spectrum we use 

(C(fc)C(fc')) = (2^)''53 (fc + fc') P{k) , P{k) = 
(((a;)2) = / dlnkAlik) ~ 0(l)A?(fc). 



This can be compared with the more precise result 
computed in [HUl EI] that uses the Green's function ap- 



This gives the result (31 ) 



fc3 



(B6) 



Alik) = AlJk) X ,m X e 



47r5 



7.5 X 10-5 



X e 



47r5 



(B13) 
(B14) 



A2(fc) ^ (C(a:)2) ^ 



a(i; • B) 



(B7) 



This expression has been plotted in figure [4] 

Now when backreaction is strong we can approximate 

E-B) 
3Hx 



/3 — 27r^Q; . , which immediately gives the approxi 



mation ( 32 1 



(B8) 



We can as well make an approximation for the case where 
/3 « 1 (negligible backreaction) and compare the result 
with the precise result (24 1, just to see how well this 



whole approximation works. For /? = 1 we have 



Alik) 



a{E ■ B) 



(B9) 



Upon using the estimate for {E ■ B) found in [501 ISB 

(BIO) 



H 

{E-B) «2.4x lO^'^^e^''^ 



and 



a : 



2H^ 



(Bll) 



we find 



Alik) 



A ff2t2 tt8 

X 5.76 X 10-8 X ^e^^^ x 



= 2.56 X 10-8 X — X — — 



2.56 X 10-8 X 



(-) 



X (27r)'* X 



4.0 X 10-5 X ,,(fc) X 



(B12) 



where in the second line we used the large ^ limit for /2. 
We infer that this quick estimate is off by a factor less 
than two. 

Actually, for some ^ the estimate comes even closer 
than this ratio ^ . Let us examine the situation at ^ = 3 
(which, for ^(7V = 60) = 2.2), corresponds to N fa 35). 
Above, we approximated the numerical function /2(C) by 
"^'^^^Y^ which yields an overestimate by a factor of 1.3. 
At the other hand, we also approximate d the numerically 
found result for {E ■ B) 



BIO 



which is an 



by the estimate 
underestimate, that for ^ = 3 only captures a fraction of 
0.73 of the true (E-B). Putting everything together one 



finds that, at ^ = 3 (N = 35), our estimate B7 with (3 set 
to one overestimates the precisely computed numerical 
result (B13l by a factor of 



1.3 



7.5 



(0.73)" 



1.3. 



(B15) 



At ^=2.2 (N=60) we find that our estimate B7 overesti- 
mates the precisely computed result by a factor of 2.5. 

Now one might introduce a fudge factor such that at 
some preferred value for ^ our approximation precisely 
matches the numerically computed result. However, we 
have just seen that the inclusion of such a fudge factor 
will induce only a small shift in our estimate that we any- 
way only trust up to corrections of order one. Besides, 
the fudge factor would always be arbitrary, as it depends 
on the preferred value of ^ where it makes both signals 
match. Therefore it seems safe to neglect it altogether. 
In Figure [7] we have for once plotted how the total power 
spectrum (including the standard slow-roll contribution) 
would shift from such a correction. In the rest of the pa- 
per we work with our uncorrected estimate for the power 
spectrum. 

N.B. This estimate involves only the gauge field contri- 
bution to the power spectrum. Apart from that there is 
always the standard slow-roll component A^ sr(^)- This 
is the dominant contribution on CMB scales. That is 
why any estimate of the total power spectrum matches 
the precise result so well on CMB scales, whatever order 
one fudge factor one chooses. 
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and therefore 




T^~2{E-Bf. (C4) 



Appendix D: Bispectrum and /jvl estimate 



The position space three point function of C, can be 
directly generaHzed from (B5l: 



10 20 30 40 50 60 



N 



FIG. 7. Evolution of the power spectrum as in Figure [4] still for 
5[Af = 60] = 2.2. The red solid line is our estimate. The blue, 
tinily dashed line is our estimate corrected with a fudge factor of 
1.3. The green, largely dashed line is our estimate corrected with 
a fudge factor of 2.5. All signals remain within an order one factor 
from the black hole bounds in dashed black. 



Appendix C: Skewness ot E ■ B 



We want to compute 



IP 

X^ 



where we used t he d efinition of the skewness (CI I 
and its estimate (C4|. {C,{xY) is positive. (Again: we 



work with negative x which yields positive {E ■ B), while 
working with x > would give {E ■ B) < 0.) 

Let us first analyze this result in the regi me w here 
backreaction is negligible, i.e. (3 = 1. Using (BIO I and 



(Bill we get 



^{(^E-B-{E-B) 



E- Bj ) -4{E ■ Bf 
E-BY), + 3(E-Bf 



(CI) 



where we used that ((^E-B^ ) ~ 2(i? -5)2 from the pre- 
vious section and in the last step we recognized that there 
are 1 + 3x2 = 7 non-connected diagrams in ( (^E ■ B^ ) , 

each one equals to (E-B)^. Using Wick's theorem we find 
many terms. All of them have the same angular integral 



d^nk,d^nk,d^nkMki)e^{k2)e*(h)e,{h)e*{k2)e*{k3) 



2^ 
3 



(C2) 



Counting all the possible pairwise contractions one finds 



3 Ojr^ 

E-B] )r = — (2/| + hhh) 
2.4-10-4^e2< 



-{E-B) 



(C3) 



(C(^) 



2— (2.4 X 10 
27 ^ 



.4.3 Hi2g6.? 



.2 X 10" 



(D2) 



Now we want to compare this with the momentum 
space bispectrum B{k), defined via 

(C(fcl)C(fc2)C(fc3)) = {27t)H^ (fci +k2+ fcs) B{ki,k2, ks) 

(D3) 

for which we can write 

3^ f d^ki f d^k2 



(27r)3 J (27r) 



B{ki,k2,-ki-k2). (D4) 



When non-Gaussianity is large mostly on equilateral tri- 
angles, the integral is supported in the region k2 ~ fci 
and 6*12 ~ 7r/3. Hence we estimate 



iaxf) = I d\ogk^k^B,^{k) 



(27r 



^fc«i?eq(fc)0(l). 



(D5) 



where B^q^k) is the bispectrum evaluated on equilateral 
triangles. Now we can compare our estimate (D2| with 
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the precisely computed result using the Green's function 
approach, that we take from result ((2.8) of [33)). 



Bcq(A:) 



1 



;(C(fcl)C(fc2)C(fc3)) 



3 X 3 X 2.8 X 10-^ iJi^gS'T? 1 



10(27r)2 ^9^6 fc6 ^ 

where we have used the large estimate 

2.8- 10^^ 



(D6) 



/3(e) 



This last result leads to 



(D7) 



(C(3 



i2n) 



gfc^Beq(fc) ~ 8.2 X 10 



TTl2 GTTi 

-12 ^ ^ I 



tHD8) 



9^6 



which agrees (surprisingly) well with ( D2 1 



In the regime of strong backreaction we can write 

{E B) Ji 1, 

/3 « — 27r^a and the estimate (|D1|) directly gives 



the generalization of ( B8 1 



(C(^) 



(D9) 



Finally we want to convert these results into a value 
for /nl- We take /^vl to be defined via 



(C(fcl)C(fe)C(fc3)) = (271)3^3 lk,+k2+h 



x(2.)^A/..A^(.)||. 

(DIO) 



This gives 



In the regime of negligible backreaction we can then take 



our estimate ( D2 1 , and conclude that 



rcq 

•I I 



2.8-10-^ e6-«A6 ^^{k) 



NL 



Ajik) 



(D13) 



This again matches the result obtained in [301 EI] by a 
more precise computation. (Of course, after that we had 
found that the expressions for {^{x)^) match so well, this 
is only a consistency check.) 

In the regime of strong backreaction, finally, we need 
to insert (D9l into ( |D12 1. Using our power spectrum 
estimate ( B8 1 we find 



_ 10(27r)2(2<)4 _ 10 
•^^^ " 9 87r2 4^3^3 " 9 ^ 



42 (D14) 



where we have used that towards the end of inflation we 
have ^ ~ 6. 



Notwithstanding the precise match between ( D2 1 and 



( D8 1 , there is a still an order one factor between the es- 
timate for the three point function (and for /nl) and 
its precisely computed numerical v alue. Again: to arrive 

for {E 



at (D2| we have used the estimate 



to arrive at ( D8 1 we have inserted 
mation 



BIO 



the 



B), and 
large ^ approxi- 



for fjii). When using precise numerical 

prescriptions rather than estimates for {E ■ B) and (f ) 
we find that our estimates overshoots the precisely com- 
puted Jml by a factor of 9.5 at ^ = 2.2 {N — 60), and 
by a factor of 3.8 at ^ = 3 (iV w 35). 

Again we will not bother introducing a fudge factor 
to close this gap at some preferred value of ^. Anyway, 
when backreaction is large /tv^ is not a suitable indicator 
for the amount of non-Gaussianity anymore. In Figure [8] 
we plot our estimate for a more meaningful quantity: the 
skewness, which is equivalent to JnlC- When backreac- 
tion becomes important, it saturates at a value of about 
one, which a posteriori justifies our approach (36 1. 



r.lO 1 1 Ilifcf Appendix E: Black hole masses 

fNL = B{k,MM^J^,^^^^, (Dll) 

In this appendix we give some details about the deriva- 
which for the equilateral case becomes tion of (III for the black hole mass and about the total 



number of efoldings enforced by a specific expansion his- 

n -B (k)"^^ 

d {ATT) L\^yK} Suppose the universe is radiation dominated right after 

(27r)^ 1 ^3 10 1 \ k^ the end of inflation. Then the expansion proceeds as 

87^2 p"^^'^^) ^ ^ "3" (27r)4 A'^(fc) 3fc3 « ~ (*Ao)^^^, so H(t) = ^. This regime starts at t^, 

^ which is not the time since the beginning of Big Bang, 

10 (27r) {Cjx) ) (D12) simply the constant = 577' where H is the Hubble 

9 87r2 A^(fc) constant at the end of inflation. We distinguish it from 
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FIG. 8. Evolution of /nl x ^ as a function of A'', for ^[A'' = 
60] = 2.2 



In other words, the black holes form after the universe 
expands by a factor since the end of inflation. The 
initial energy stored inside the volume H^^e^ was Af^y — 
10 e"^^ g, but during this extra expansion it scales down 
(redshifts) by the factor e^, so it becomes 

Mbh^ lOe^^g. 

It should be stressed that specifying the energy density at 
the end of inflation and at reheating, directly determines 
the number of efoldings corresponding to any scale (and 
in particular CMB scales) according to [51] 



N{k) = 62 - log 



k 



log 



lOi^GeV 



1/4 



(El) 



■log 



1/4 



1/4 



1 V'^' 

3 ^ 

Preh 



-1„N 



the decreasing H{t) = The wavelength lt„ = H 
grows as k = H~^{t/tof/'^e'^ = H-^{2Ht)^^^e'^ . The 
horizon size 1/H(t) = 2t grows and becomes equal to It 
(and black holes form) at 



i.e. at 



2t ■ 



(2Ht) 



iJ-i(2i7i)i/2e^ 



1/2 _ 



(t/to) 



1/2 ^ ^N^ 



where is the energy density during inflation when the 
mode k left the horizon. Vend is the energy density at the 
end of inflation, preh is the energy density at reheating 
and the subscript refers to todays value. Taking for 
example p^eh = Vend = and Vk = m^l5'^M^/2 

with m = 6 X lO^Mp gives Ncmb = N{aoHo) ~ 64. We 
use this value in our discussion of primordial black holes, 
but since the difference between 60 and 64 changes very 
little in our numerics, for simplicity we use Ncmb — 60 
in the rest of the paper. 



[1] R. Kallosh and A. Linde, "New models of chaotic 
inflation in supergravity," JCAP 1011, Oil (2010) 
[arXiv: 1008.3375 [hep-th]]. 

[2] R. Kallosh, A. Linde and T. Rube, "General inflaton po- 
tentials in supergravity," Phys. Rev. D 83, 043507 (2011) 
[arXiv: 101 1.5945 [hep-th]]. 

]3] R. Kallosh, A. Linde, K. A. Olive and T. Rube, "Chaotic 
inflation and supersymmetry breaking," Phys. Rev. D 84, 
083519 (2011) ]arXiv:1106.6025 [hep-th]]. 

j4] M. Kawasaki, M. Yamaguchi and T. Yanagida, "Natural 
chaotic inflation in supergravity," Phys. Rev. Lett. 85, 
3572 (2000) [arXiv:hep-ph/0004243[. 

[5] M. Yamaguchi, "Natural double inflation in super- 
gravity," Phys. Rev. D 64, 063502 (2001) ]arXiv:hep- 
ph/0103045]. 

}6] M. Kawasaki and M. Yamaguchi, "A supersymmetric 
topological inflation model," Phys. Rev. D 65, 103518 
(2002) ]arXiv:hep-ph/0112093]. 

]7] M. Yamaguchi and J. Yokoyama, "Chaotic hybrid new 
inflation in supergravity with a running spectral index," 
Phys. Rev. D 68 (2003) 123520 [arXiv;hep-ph/0307373[. 

[8[ P. Brax and J. Martin, "Shift symmetry and infla- 
tion in supergravity," Phys. Rev. D 72, 023518 (2005) 
[arXiv:hep-th/0504168]. 

[9] R. Kallosh, "On inflation in string theory," arXiv:hep- 
th/0702059. 

[10[ K. Kadota and M. Yamaguchi, "D-term chaotic infla- 



tion in supergravity" Phys. Rev. D 76, 103522 (2007) 
[arXiv:0706.2676 [hep-ph[[. 

[11[ S. C. Davis and M. Postma, "SUGRA chaotic infla- 
tion and moduli stabihsation," JCAP 0803, 015 (2008) 
[arXiv:0801.4696 [hep-ph[[. 

[12[ M. B. Einhorn and D. R. T. Jones, "Inflation with Non- 
minimal Gravitational Couplings in Supergravity," JHEP 
1003, 026 (2010) [arXiv:0912.2718 [hep-ph[[. 

[13[ S. Ferrara, R. Kallosh, A. Linde, A. Marrani and 
A. Van Proeyen, "Jordan Frame Supergravity and In- 
flation in NMSSM," Phys. Rev. D82, 045003 (2010) 
[arXiv: 1004.0712 [hep-th]|. 

[14[ H. M. Lee, "Chaotic inflation in Jordan frame supergrav- 
ity," JCAP 1008, 003 (2010) [arXiv: 1005.2735 [hep-ph[[. 

[15[ S. Ferrara, R. Kallosh, A. Linde, A. Marrani and A. Van 
Proeyen, "Superconformal Symmetry, NMSSM, and In- 
flation," arXiv: 1008.2942 [hep-th[. 

[16[ F. Takahashi, "Linear Inflation from Running Kinetic 
Term in Supergravity," arXiv: 1006.2801 [hep-ph[; 

[17[ K. Nakayama and F. Takahashi, "Running Kinetic Infla- 
tion," arXiv:1008.2956 [hep-ph[. 

[18[ E. Silverstein and A. Westphal, "Monodromy in the 
CMB: Gravity Waves and String Inflation," Phys. Rev. 
D 78, 106003 (2008) [arXiv:0803.3085 [hep-th[[. 

[19[ L. McAllister, E. Silverstein and A. Westphal, "Grav- 
ity Waves and Linear Inflation from Axion Monodromy," 
arXiv:0808.0706 [hep-th[. 



17 



[20] M. Yamaguchi and J. Yokoyama, "New inflation in su- 
pergravity with a chaotic initial condition," Pliys. Rev. 
D 63, 043506 (2001) [arXiv:hcp-ph/0007021]. 

[21] R. Kallosh and A. D. Linda, "Landscape, the scale of 
SUSY breaking, and inflation," JHEP 0412, 004 (2004) 
[arXiv:hep-th/0411011]; J. J. Blanco-Pillado, R. Kallosh 
and A. D. Linde, "Supersymmetry and stability of 
flux vacua," JHEP 0605, 053 (2006) [hep-th/0511042]; 
R. Kallosh and A. D. Linde, "Testing String Theory with 
CMB," JCAP 0704, 017 (2007) [arXiv:0704.0647 [hep- 
th]]. 

}22] A. Linde, Y. Mambrini and K. A. Olive, "Supersymmetry 

Breaking due to Moduli Stabilization in String Theory," 
Phys. Rev. D 85, 066005 (2012) ]arXiv: 111 1.1465 [hep- 
th]]. 

}23] E. Dudas, A. Linde, Y. Mambrini, A. Mustafayev and 
K. A. Olive, "Strong moduli stabilization and phe- 
nomenology," arXiv: 1209.0499 [hep-ph]. 

]24] S. Ferrara and R. Kallosh, "Creation of Matter in the Uni- 
verse and Groups of Type E7," JHEP 1112, 096 (2011) 
JarXiv: 11 10.4048 ]hep-th]]. 

125] A. D. Linde and V. Mukhanov, "Nongaussian isocur- 
vature perturbations from inflation," Phys. Rev. D 
56, 535 (1997) [arXiv:astro-ph/9610219]; K. Enqvist 
and M. S. Sloth, "Adiabatic CMB perturbations in 
pre big bang string cosmology," Nucl. Phys. B 626, 
395 (2002) [arXiv:hcp-ph/0109214]; D. H. Lyth and 
D. Wands, "Generating the curvature perturbation with- 
out an inflaton," Phys. Lett. B 524, 5 (2002) ]arXiv:hcp- 
ph/0110002]; T. Moroi and T. Takahashi, "Effects of 
cosmological moduli fields on cosmic microwave back- 
ground," Phys. Lett. B 522, 215 (2001) ]Erratum-ibid. 
B 539, 303 (2002)] ]arXiv:hcp-ph/0110096]. 

126] V. Dcmozzi, A. Linde and V. Mukhanov, "Supercurva- 
ton," JCAP 1104, 013 (2011) JarXiv: 1012.0549 Jhep-th]]. 

127] M. M. Anber and L. Sorbo, "N-flationary magnetic 
fields," JCAP 0610, 018 (2006) lastro-ph/0606534]. 

128] M. M. Anber and L. Sorbo, "Naturally inflating on steep 
potentials through electromagnetic dissipation," Phys. 
Rev. D 81, 043534 (2010) larXiv:0908.4089 [hep-th]]. 

129] R. Durrer, L. Hollenstein and R. K. Jain, "Can slow roll 
inflation induce relevant helical magnetic fields?," JCAP 
1103, 037 (2011) [arXiv:1005.5322 lastro-ph.CO]]. 

130] N. Barnaby and M. Pcloso, "Large Nongaussianity in 
Axion Inflation," Phys. Rev. Lett. 106 (2011) 181301 
larXiv: 101 1.1500 lliep-pli]]. 

131] N. Barnaby, R. Namba and M. Peloso, "Phenomenology 
of a Pseudo-Scalar Inflaton: Naturally Large Nongaus- 
sianity," JCAP 1104, 009 (2011) larXiv: 1102.4333 lastro- 
ph.CO]]. 

132] N. Barnaby, E. Pajcr and M. Pcloso, "Gauge Field 
Production in Axion Inflation: Consequences for Mon- 
odromy, non-Gaussianity in the CMB, and Gravitational 
Waves at Interferometers," Phys. Rev. D 85, 023525 
(2012) [arXiv:1110.3327 lastro-ph.CO]]. 

133] P. D. Meerburg and E. Pajer, "Observational Con- 
straints on Gauge Field Production in Axion Inflation," 
arXiv: 1203.6076 [astro-ph.CO]. 

134] J. L. Cook and L. Sorbo, "Particle Production during In- 
flation and Gravitational Waves Detectable by Ground- 
Based Interferometers," Phys. Rev. D 85 (2012) 023534 
lErratum-ibid. D 86 (2012) 069901] larXiv:1109.0022 
[astro-ph.CO]]. 

135] M. M. Anber and L. Sorbo, "Non-Gaussianities and 



Chiral Gravitational Waves in Natural Steep Inflation," 
Phys. Rev. D 85 (2012) 123537 larXiv:1203.5849 lastro- 
ph.CO]]. N. Barnaby, J. Moxon, R. Namba, M. Pcloso, 
G. Shiu and P. Zhou, "Gravity Waves and Non-Gaussian 
Features from Particle Production in a Sector Gravita- 
tionally Coupled to the Inflaton," Phys. Rev. D 86 (2012) 
103508 larXiv:1206.6117 lastro-ph.CO]]. 

[36] A. S. Josan, A. M. Green and K. A. Mahk, "Gener- 
alised constraints on the curvature perturbation from 
primordial black holes," Phys. Rev. D 79 (2009) 103520 
[arXiv:0903.3184 [astro-ph.CO[]. 

[37] B. J. Carr, K. Kohri, Y. Sendouda and J. 'i. Yokoyama, 
"New cosmological constraints on primordial black holes," 
Phys, Rev. D 81, 104019 (2010) larXiv:0912.5297 lastro- 
ph.CO]]. 

[38] C. T. Byrnes, E. J. Copeland and A. M. Green, 
"Primordial black holes as a tool for constraining 
non-Gaussianity," Phys. Rev. D 86, 043512 (2012) 
[arXiv:1206.4188 [astro-ph.CO|]. 

[39] E. Bugaev and P. Klimai, "Curvature perturbation 
spectra from waterfall transition, black hole con- 
straints and non-Gaussianity," JCAP 1111, 028 (2011) 
larXiv: 1107.3754 lastro-ph.CO]]. E. Bugaev and P. Kli- 
mai, "Formation of primordial black holes from non- 
Gaussian perturbations produced in a waterfall transi- 
tion," Phys. Rev. D 85, 103504 (2012) larXiv:1112.5601 
lastro-ph.CO]]. P. A. Klimai and E. V. Bugaev, "Primor- 
dial black hole formation from non-Gaussian curvature 
perturbations," arXiv:1210.3262 [astro-ph.CO]. 

[40] D. H. Lyth, "The hybrid inflation waterfall and the pri- 
mordial curvature perturbation," JCAP 1205, 022 (2012) 
larXiv:1201.4312 lastro-ph.CO]]. 

[41] S. Shandera, A. L. Erickcek, P. Scott and J. Y. Galarza, 
"Number Counts and Non-Gaussianity," arXiv:1211.7361 
[astro-ph.CO]. 

[42] P. Fayet, "Phenomenology Of Supersymmetry," Talk at 
the XVIIth Rencontre de Moriond, Ecole Normale Su- 
perieure preprint LPTENS 82/10 (1982); S. Weinberg, 
"Cosmological Constraints on the Scale of Supersymme- 
try Breaking," Phys. Rev. Lett. 48, 1303 (1982); J. R. El- 
lis, A. D. Linde and D. V. Nanopoulos, "Inflation Can 
Save the Gravitino," Phys. Lett. B 118, 59 (1982); 
L. M. Krauss, "New Constraints on Ino Masses from Cos- 
mology. 1. Supersymmetric Inos," Nucl. Plij's. B227, 556 
(1983); J. R. Elhs, J. S. Hagchn, D. V. Nanopoulos, 
K. A. Olive and M. Srcdnicki, "Supersymmetric Relics 
from the Big Bang," Nucl. Phys. B 238, 453 (1984); 
M. Y. Khlopov, A. D. Linde, "Is It Easy to Save the 
Gravitino?," Phys. Lett. B138, 265-268 (1984). 

[43] E. Komatsu et al. IWMAP Collaboration], "Seven- Year 
Wilkinson Microwave Anisotropy Probe (WMAP) Ob- 
servations: Cosmological Interpretation," Astrophys. J. 
Suppl. 192, 18 (2011) larXiv:1001.4538 [astro-ph.CO]]. 

[44] T. Banks, M. Dine, P. J. Fox and E. Gorbatov, "On 
the Possibility of Large Axion Decay Constants," JCAP 
0306 (2003) 001 [hep-th/0303252]. 

[45] C. -M. Lin and K. -W. Ng, "Primordial Black Holes from 
Passive Density Fluctuations," arXiv: 1206. 1685 [hcp-pli[. 

[46] A. D. Linde, "Chaotic Inflation," Phys. Lett. B 129, 177 
(1983). 

[47] A. D. Linde and V. Mukhanov, "The Curvaton Web," 
JCAP 0604, 009 (2006) larXiv:astro-ph/0511736]. 

148] A. D. Linde, "Eternally Existing Self-reproducing 
Chaotic Inflationary Universe," Phys. Lett. B 175, 395 



18 



(1986). 

[49] V. F. Mukhanov and G. V. Chibisov, "Quantum Fluctu- 
ation And Nonsingular Universe," JETP Lett. 33, 532 

(1981) [Pisma Zh. Eksp. Teor. Fiz. 33, 549 (1981)]; 
S. W. Hawking, "The Development Of Irregularities In A 
Single Bubble Inflationary Universe," Phys. Lett. B 115, 
295 (1982); A. A. Starobinsky, "Dynamics Of Phase Tran- 
sition In The New Inflationary Universe Scenario And 
Generation Of Perturbations," Phys. Lett. B 117, 175 

(1982) ; A. H. Guth and S. Y. Pi, "Fluctuations In The 
New Inflationary Universe," Phys. Rev. Lett. 49, 1110 
(1982); J. M. Bardeen, P. J. Steinhardt and M. S. Turner, 
"Spontaneous Creation Of Almost Scale - Free Density 
Perturbations In An Inflationary Universe," Phys. Rev. 
D 28, 679 (1983); V. F. Mukhanov, "Gravitational In- 
stability Of The Universe Filled With A Scalar Field," 
JETP Lett. 41, 493 (1985) ]Pisma Zh. Eksp. Teor. Fiz. 
41, 402 (1985)]; M. Sasaki, "Large Scale Quantum Fluc- 
tuations in the Inflationary Universe," Prog. Theor. Phys. 



76, 1036 (1986)." 
[50] N. Arkani-Hamed and S. Dirnopoulos, "Supersymmet- 
ric unification without low energy supcrsymmctry and 
signatures for fine-tuning at the LHC," JHEP 0506, 
073 (2005) [arXiv:hep-th/0405159]; G. F. Giudice and 
A. Strumia, "Probing High-Scale and Split Supersymme- 
try with Higgs Mass Measurements," Nucl. Phys. B 858, 
63 (2012) ]arXiv:1108.6077 ]hep-ph]]; M. Ibc, S. Mat- 
sumoto and T. T. Yanagida, "Pure Gravity Mediation 
with w-a/a = 10-lOOTeV," Phys. Rev. D 85, 095011 (2012) 
JarXiv: 1202.2253 Jhep-ph]]; E. Dudas, A. Linde, Y. Mam- 
brini, A. Mustafayev and K. A. Olive, "Strong moduli 
stabilization and phenomenology," arXiv:1209.0499 [hep- 
ph]. 

[51] A. R. Liddle and D. H. Lyth, "The Cold Dark Matter 
Density Perturbation," Phys. Rept. 231 (1993) 1 [astro- 
ph/9303019[. 



